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Abstract 



We introduce an exactly solvable model to study decoherence of a central spin 
interacting with a spin bath where the coupling is mediated by phonons which we 

■ assume to be in a coherent state or thermal distribution. For the coherent state case, 
we find that the decoherence factor decays in a Gaussian fashion and it becomes 
independent of the phonon frequencies at short times. If the phonon energies are 

■ much larger than spin-phonon coupling or bath spins are fully polarized, decoher- 
^ |. ence time becomes independent of the initial phonon state. For the thermal state 

case, phonons play more important role in decoherence with increasing tempera- 
ture. We also discuss possible effects of the temperature on spin bath contribution 
to decoherence. 
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Decoherence is the key concept for understanding the emergence of classical 
states out of a quantum system [I] . This phenomenon is also the main challenge 
^ I in quantum information processing [2]. Coupling of central two level system 

o3 I (qubit) to environment, leads to loss of phase relation between the states 

of the qubit. Therefore, the superposition of the qubit states evolves into a 
statistical mixture of the states, so called pointer states. These states are 
determined by the form of the system-environment interaction. If the qubit 
starts from a pure state where it is decoupled from the environment, in time, 
the qubit and the environment become quantum mechanically correlated. As 
the qubit gets entangled with the environment, it can no longer be described 
by a pure state. Although decoherence concept seems to solve most of the 
puzzle of the emergence of classical! ty, there are still open questions like "How 
does the information flow from system to environment?" . Understanding this 
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information transfer is believed to be crucial for explaining the objectivity of 
the classical world [3P] . 



Being the elementary quantum information units, qubits are one of the most 
extensively studied open quantum systems. Especially solid state qubits (quan- 
tum dots, SQUIDs, magnetic molecules, etc.) have attracted a great interest 
due to their scalability which is an indispensable criterion for realistic quantum 
information processing. However, most important drawback of these systems is 
their relatively strong couplings to the environment. Understanding the mech- 
anism of these interactions is crucial for the implementation of error correction 
techniques [5] and/or error avoiding strategies [6]. Starting with the pioneer- 
ing works of Caldeira and Leggett [7118] . the crucial effects of environment on 
the dynamics of the central system has been studied with different models. 
Among them spin-boson model has attracted much attention P|ITU] . Now, it is 
a well understood environmental model. However, this model is inadequate in 
most situations where localized environmental modes act as a main source of 
decoherence [TT][T2] . In these cases spin bath models are used to describe the 
environment. In spite of numerous theoretical works, including both analyt- 
ical approaches P^14|15|16lll7j and numerical simulations p!5|19ll2U|[^ . spin 
bath decoherence is still a hot subject. This is due to the rich dynamics of 
spin models with different intra-bath couplings. Quantum dots are extensively 
studied systems, both theoretically and experimentally, where the hyperfine 
interaction with nuclear spins is dominating mechanism of decoherence |22] . 

Generally, bosonic and fermionic modes are considered to be effective at dif- 
ferent time scales and they are coupled to the qubit independently. However, 
this is not always the case. For instance, recent experiments on particular sin- 
gle molecular magnets show that these two mechanisms cooperate together 
[22121] ■ It has been proposed that the Waller mechanism, modulation of the 
dipolar fields by atomic vibrations, can play an important role [23|. Phonon 
assisted hyperfine interaction in quantum dots is another example [2S]- Devi- 
ations of the nucleus positions due to lattice vibrations modify the hyperfine 
coupling with electron spin. A final example can be given from optical lat- 
tices where coupling strengths between spins trapped deep inside a confining 
potential change with lattice oscillations. Starting with the theoretical study 
of Jaksch et al. [27J, ultra-cold atomic gases in optical lattices have attracted 
great attention. Possibility of controlling the interactions among trapped par- 
ticles is most advantageous property of these systems. This peculiarity enables 
to mimic various spin models such as Ising, XY, Heisenberg and so on (see 
for review [2B])- These developments lead to various applications in quantum 
information processing [21] and study of spin bath decoherence in a controlled 
way [30|31f32j . 




Inspired by these observations, we introduce a pure dephasing model where 
the interaction of the central two-level system with environmental spins is 
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mediated by phonons. Study of pure dephasing model is motivated by two 
observations. Firstly, dissipative processes where the energy exchange occurs 
between subsystems have typically longer time scales than pure dephasing 
processes Secondly, exact solubility of the model gives a more clear un- 
derstanding of the decoherence process. We neglect the self-Hamiltonians of 
the central spin and the spin bath. In particular, we don't consider any in- 
teraction among the bath spins. We assume that low energy physics of the 
system is governed by the effective Hamiltonian 

N N 

H = C^Y1 [^Ok + ^k [pI + Pk) Skz + ^kPlPk (1) 
k=l k=l 

where Cz and Skz are z-components of the Pauli spin operators for central two- 
level system and k^^ spin of the bath, respectively. is the total number of 
environmental spins. We are using units such that the Planck and the Boltz- 
mann constants are unity, pi and pk are the boson creation and annihilation 
operators with commutation relation [pfc,p|./] = Sk,k'- Energy eigenvalues of the 
phonon bath are Qk, and the coupling strengths are uok and Uk- Our model 
is similar to the one proposed by Zurek where the central two level system 
is directly coupled to spin bath In our model this coupling occurs with 
the help of oscillatory modes. When Uk and Qk vanish our model reduces to 
Zurek's. The model Hamiltonian describes a system where the interaction be- 
tween the central two-level system and the bath spins is distance dependent 
and this distance is modified by some vibrational modes. 

First, we solve the case where the qubit is surrounded by spins almost lo- 
calized at different positions, for example at lattice points of a solid. The 
interaction strengths between the system and a bath spins change with the 
distance between them. Considering the displacement of these atomic posi- 
tions as macroscopic vibrations, we model them by coherent states which are 
the most classical states of phonons. An atom, confined in a harmonic po- 
tential, satisfies the minimum position-momentum uncertainty when it is in a 
coherent state which is nothing but a Gaussian wave function displaced from 
the origin. Furthermore, it oscillates while preserving its shape, i.e. it remains 
as a coherent state. We assume that initially the system and the environment 
are uncorrelated so that the initial wave function can be written as a product 
state, 

N 

1^(0)) = (ctI T) + q| i)) (g) U) + Pk\ Ik)) \\k) (2) 

k=l 

where | t)(| ]k)) and | |)(| 1^)) are normalized eigenstates of Cz{skz) with 
eigenvalues +1 and-1, respectively. Expansion coefficients satisfy |c|p + |cjp = 
|«fcP + lAP = 1 so that l^'(O)) is normalized. \Xk) is the coherent state cor- 
responding to the annihilation operator pk with eigenvalue \k so that Pk\^k) = 
Ajt|A/c). With the help of the harmonic displacement operators D{a) = e"^^""*^, 
Hamiltonian can be diagonalized easily. Applying the propagator e~**^, we 
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can calculate the time evolution of the wave function which can be written as 
\^{t)) = ctI t)|i^+W) + q| i)\B4t)) where 



N 

k=l 



B±{t)) = (g) {a^A^l U)H)+PkAt\ h)\ut)) . (3) 



Here \u^) are the coherent states with eigenvalues 
and 

^± ^g*s\V ^^fc J ^=pjf<^oj^^T«s^(Re[Afc]sin{nfct)+Im[Afc](l-cos(nfct))) 

Total density matrix is given by p = Reduced density matrix of 

the central system is obtained by tracing over all the environmental degrees 
of freedom as Pc = TrbathP- In c^-basis, the reduced density matrix is given by 



Pc 



* * I 1 2 

|q| j 



(6) 



Magnitude of the off-diagonal matrix element is determined by the decoherence 
factor 

TV 

r{t) = UMA^'^ti^klO + WAfA^iutlui)) (7) 
fc=i 

which can be written more explicitly as 



^ -Q ^-4HT(l-cos(nfct)) ^|^^|2^_i2a;ofct-i4^(Re[Afc]sin(nfct)+Im[Afc](l-cos(f2fct))) 
fe=l 



_|_j^^|2g«2a)0fe«+»4j^(Re[Afe]sin(ai;t)+Im[Afc](l-cos(Ofct)))^ 

(8) 



At t — 0, r — 1 and as t increases, in general, it decays to zero which means 
that interference of the states | f) and | i) is totally suppressed. At short 
enough times we can expand the trigonometric functions by treating Qkt^s as 
small parameters to obtain 

TV 

~ Yl e"^'^**^ (^|Q,^|2g-«t{4a;fcRe[Afc]+2a;ofc) _^ |^^|2gjt(4u;fcRe[Afc]+2a;ofe)^_ ^g^ 
fc=l 

If either the coupling strengths a;fc's and a;ofc's or coherent state eigenvalues 
Afc's are random enough, the second factor in the product leads to further 
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suppression of the coherence factor so that r decays in Gaussian form for 
large N [35] . 



r{t)\ ^ e 



1 2 ^ J8 1 afe 1 2 1 & 1 2 (2a;fc Re [ Afc ] +u;„fe ) 2 +2a;2 ) 



(10) 



Therefore phonon energies do not play any role for short time decoherence 
of the central system. The decoherence time is determined by the coupling 
constants and the initial configurations of the phonons and spin bath states. 
It is interesting to note that even if all the bath spins are polarized in one 
direction, i.e. |afcp = 1, system still losses its coherence. This behavior is a 
result of presence of the phonons in the environment. It is also interesting 
that phonon state eigenvalues (A^'s) do not affect the decoherence time in this 
case. It is obvious that in the limit of fi^ — > and uj^ 0, our Hamiltonian 
is reduced to Zurek's model where decoherence is due to direct spin-spin in- 
teractions only without phonon contribution. In this case initial configuration 
of the spin bath becomes crucial. 

Another interesting case is the VLk/oJk oo limit where the decoherence factor 
r{t) = nLi(l"fePe"*^'^o'=* + |/5fcpe*2'^o'=*). Therefore, r{t) depends on the initial 
configuration of bath spins only and it becomes independent of the initial 
phonon state eigenvalues. Since the separation of energy levels of phonons be- 
comes very high, phonon states do not change in time and remain uncorrelated 
to system and bath spins. 

Now, we analyze the case where phonons are in a thermal equilibrium rather 
than a coherent state. Such a situation can physically be realized when the 
atoms carrying bath spins are brought in contact with a heath bath to ther- 
malize before t = 0. For thermal states phonon density matrix is given by 



We assume that the bath spins are in a separable state at t = as before. Since 
in the Hamiltonian there are no intra-bath terms for the spins, heath bath 
thermalizing the phonons will simply randomize the initial spin directions. 
As we shall discuss below, if bath spins have individual energy levels for up 
and down configurations, heath bath will determine the initial occupation 
numbers for the two possible states in accordance with the Gibbs factors. 
Time evolution of the total density matrix is given by p{t) = e~*^*p(0)e*''^*. 
Using the over-completeness relation 



Pp(0) = (g)(l-e-^) 




nk){nk 




k nfc=0 
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and the number state representation of coherent states 



{n\X)=e-\^\''^^. (13) 



it is straight forward to calculate the reduced density matrix of the central 
system. In this case decoherence factor becomes 



2 

^ -4^(l-cos(nfci))coth(^) 



r(t) = []e ^ ^ (|afc|2e-'2"«'=* + (14) 



k=l 



We first note that for Qk/T —>■ oo, Eq. f|T^ and Eq. ([8]) become identical pro- 
vided that Afc = 0. This is a consistency check for two phonon states, coherent 
states and thermal states, that we have discussed because at low temperatures 
thermal state approaches the ground state of the harmonic oscillators which 
are nothing but the coherent states with vanishing eigenvalues. 

According to Eq. ([HD, decoherence factor has two contributions, coming from 
phonons and spins. The two mechanisms act simultaneously in decoherence 
of the central spin. Depending upon the interaction strengths, one of them 
can become the dominant mechanism. At very low temperatures, where the 
hyperbolic cotangent term is approximately unity, the first term becomes in- 
dependent of flk values provided that t is small enough. For large tempera- 
tures, decoherence factor becomes an exponentially decaying function of T. 
It is possible to generalize the model Hamiltonian by adding a s^^-dependent 
intra-bath term for individual spins. In this case the heath bath will not only 
thermalize the phonons but also it will determine the |afcP/|/?fcP ratio. For 
example, at very large temperatures the ratio will tend to unity and hence the 
spin bath will have a more important contribution to decoherence in compar- 
ison to lower temperatures. 

In conclusion, we examined a spin decoherence model where the interaction 
with the bath spins are modified by phonons. Coherent states of phonons 
correspond to almost localized bath spins. In this case we find that initial 
decoherence rate does not depend on the phonon energies. Furthermore, for 
polarized bath spins it becomes independent of the initial phonon states. Ther- 
mal phonon distribution is the other case where we found an explicit solution 
of the model Hamiltonian. At high temperatures phonons play a more impor- 
tant role in the dephasing process. 
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